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Abstract 

The density dependence of the symmetry energy, instrumental in understanding the behaviour of the 
asymmetric nuclear matter, is investigated within the extended relativistic mean field (ERMF) model which 
includes the contributions from the self and mixed interaction terms for the scalar-isoscalar (cr), vector- 
isoscalar {a>) and vector-isovector (p) mesons upto the quartic order Each of the 26 different parameter- 
izations of the ERMF model employed are compatible with the bulk properties of the finite nuclei. The 
behaviour of the symmetry energy for several parameter sets are found to be consistent with the empirical 
constraints on them as extracted from the analyses of the isospin diffusion data. The neutron-skin thickness 
in the ^''^Pb nucleus for these parameter sets of the ERMF model lie in the range of ~ 0.20 - 0.24 fm which 
is in harmony with the ones predicted by the Skyrme Hartree-Fock model. We also investigate the role of 
various mixed interaction terms which are crucial for the density dependence of the symmetry energy. 
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I. INTRODUCTION 



The accurate knowledge of the equation of state (EOS) for asymmetric nuclear matter is im- 
portant in understanding the structure of finite nuclei away from the stability line and critical 
issues in astrophysics. The EOS for asymmetric nuclear matter is mainly governed by the density 
dependence of the nuclear symmetry energy. Indeed, significant progress has been made in under- 
standing the behaviour of the symmetry energy at the subnormal densities from the analyses of the 
isospin diffusion data in heavy ion collision l^^M and from the available data for the neutron skin 
thickness of several nuclei [j5J. The constraints on the density dependence of the nuclear symmetry 
energy as extracted from the isospin diffusion data agree with the ones deduced from the isoscaling 
analyses of isotope ratios in intermediate energy heavy-ion collisions [16]. The experimental data 
on the isotopic dependence of the nuclear giant monopole resonance in even-A Sn isotopes 13, [80 
also provides some informations on the nuclear symmetry energy which is in agreement with those 
derived from the analyses of the isospin diffusion data. The behaviour of nuclear symmetry energy 
at supranormal densities is largely unknown. Theoretically, at supranormal densities, even the is- 
sue of whether the symmetry energy increases or decreases with density still remains unresolved. 
The precise measurements of the properties of the compact stars and the transport model analyses 
of the heavy-ion collisions at intermediate and high energies can provide some constraints on the 
high density behaviour of the symmetry energy. 

The density dependence of symmetry energy obtained using the Skyrme Hartree-Fock models 
have been confronted with the constraints extracted by analyzing the isospin diffusion data Is]. 
The symmetry energy and its slope and the curvature parameter at the saturation density for only 
four out of 21 different parameter sets of the Skyrme interactions are found to be consistent with 
the ones extracted from the isospin diffusion data [a] . It may be pointed out that all of these four 
parameter sets are obtained by fitting the experimental data for the binding energies and charge 
radii for finite nuclei. Scenario is not the same when the similar investigation nl(M is carried out 
using three different versions of the relativistic mean field (RMF) models, namely, (i) models with 
meson field self interaction (ii) models with density dependent meson-nucleon couplings and (iii) 
point coupling models without meson fields. Out of 23 parameter sets of these RMF models, only 
a few are found to yield symmetry energies and their density dependence which are consistent 
with the empirical constraints imposed by isospin diffusion data. In particular, 10 different pa- 
rameter sets were considered for the model (i), but, only two of them could yield behaviour of 
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the symmetry energy consistent with the empirical constraints. Both of these parameter sets are 
obtained using nuclear matter observables, instead of fit to the bulk properties of finite nuclei. The 
density dependence of the symmetry energy is studied also for several other parameter sets of the 
RMF model [11]. None of these parameter sets yield acceptable results for the density depen- 
dence of the symmetry energy. Very recently [|l2n . density dependence of the symmetry energy 
has been studied using extended relativistic mean-field (ERMF) model which includes the contri- 
butions from self and mixed interaction terms for the cr, co and p mesons upto the quartic order. 
Even in the ERMF model, the density dependence of the symmetry energy obtained for several 
para meter sets are found to be inadequate. One of the existing parameter set is fine tuned in Ref. 
II12I1 so that the resulting behaviour of the symmetry energy can fulfill the empirical constraints on 
them. However, this parameter set is not capable of reproducing the experimental data on the bulk 
properties of finite nuclei. We would like to emphasize that the ERMF model, due to the presence 
of the various mixed interaction terms, can yield wide variations in the density dependence of 
the symmetry energy without affecting the quality of the fit to the bulk properties of finite nuclei 
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l4 \. Nevertheless, the slope of the symmetry energy at the saturation density or alternatively 



the neutron-skin thickness for the parameter sets of the ERMF model considered in Ref. [|l2l] are 
either too low or quite high. 

In the present work we investigate the density dependence of the symmetry energy using 26 
different parameterizations of the ERMF model. All the parameterizations of the ERMF model 
considered are obtained by fitting the experimental data on the binding energy and charge radius 
for the finite nuclei. Furthermore, these parameter sets yield the neutron-skin thickness in ^ospj-, 
nucleus which vary over a wide range from 0.16 - 0.28 fm. We find that quite a few of these 
parameterizations can fulfill the empirical constraints on the symmetry energy. Role of various 
mixed interaction terms are also investigated. 

The paper is organized as follows. In Sec. II we describe the ERMF model in brief. In Sec. 
Ill, we provide the expressions used to compute various quantities associated with the nuclear 
matter along with the empirical constraints on them. In Sec. IV, the results obtained using dif- 
ferent parameterizations of the ERMF model are confronted with the empirical constraints on the 
symmetry energy as extracted from the analyses of the isospin diffusion data. The role of mixed 
interaction terms of the ERMF model which are important in determining the variations in the 
density dependence of the symmetry energy is investigated in Sec. V. In Sec. VI we state our 
conclusions. 
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II. EXTENDED RELATIVISTIC MEAN FIELD MODEL 



The ERMF model includes the contributions from the self and mixed interaction terms for the 
scalar-isoscalar (cr), vector-isoscalar (co) and vector-isovector (p) mesons upto the quartic order. 
Mixed interaction terms involving p-meson field enables one to vary the density dependence of 
the symmetry energy coefficient and the neutron skin thickness in heavy nuclei over a wide range 



1411 . The contribution from the self 



without affecting the other properties of finite nuclei [|13L 
interaction of a»-mesons plays important role in varying the high density behaviour of the EOS and 



also prevents instabilities in the calculation of the EOS IIS 



1611 . On the other hand expectation 



value of the p-meson field is order of magnitude smaller than that for the a»-meson field [|l7ll . Thus, 
inclusion of the p-meson self interaction can affect the properties of the finite nuclei and neutron 

n 

stars only very marginally [11611 . The effective Lagrangian density for the ERMF model can be 
written as, 

-C = -CnM + -Ccr + -Cio + -Cp + -Ca-ajp- (1) 

where the nucleonic and mesonic Lagrangian £.nm can be written as. 



(2) 



J=n,p 



Here, the sum is taken over the neutrons and protons, r are the isospin matrices. The Lagrangian 
describing self interactions for cr, co, and p mesons can be written as. 



2 2 4 

m^cr , 



1 



1 



-Cp = -^rPpyp^" + ::m P^ff. 



(3) 

(4) 
(5) 



The ct/"', p'"' are field tensors corresponding to the io and p mesons, and can be defined as a/"' = 
&'oj'' - d'uf and (f^ = d^p" - d^p^. The mixed interactions of cr, co, and p mesons Xo-^p can be 
written as. 



V\ ^ ...2 ,u, V2 2 2 2,, ,u, „ 2 



'^^g>Wp,p-. 



+ -^gim'cr^p^p^ + 



(6) 



4^2 "5 0- i-f"- AM^ ^ 
The JLem is Lagrangian for electromagnetic interactions and can be expressed as. 



(7) 



where, A is the photon filed and F'^'^ = &^A^ - B'A^. The equation of motion for nucleons, mesons 
and photons can be derived from the Lagrangian density defined in Eq.©- The contributions from 
Eq. ([7]) are included only for the case of finite nuclei. 



III. EMPIRICAL CONSTRAINTS ON SYMMETRY ENERGY 



The symmetry energy ^sym, slope L and curvature ^sym can be evaluated as 



^sym(p) ~ ^ 



1 SE{p,S) 



L = 3p( 



2 J(52 
dE,yM 



6=0 



sym 



9pI 



dp 



dp^ 



P=PO 



P=PO 



(8) 

(9) 

(10) 
(11) 



where, po is the saturation density, E(p, S) is the energy per nucleon at a given density p and 
asymmetry 6 = {p„ -pp)/p. The density dependence of the symmetry energy can also be expressed 
in terms of E^yjip), L and ^sym as, 

,2 



Esymip) - ■E'sym(Po) + L 



P-Po 



+ 



sym Ip — Po 



(12) 



3po / 2 \ 3p() 

The above equation represent very well the behaviour of the symmetry energy at subnormal den- 
sities. At supranormal densities, one needs to include in Eq. (fT2)) the contributions of higher order 
terms lllSn . We also evaluate, 



Ko = 9pl 



d^Eoip) 

Jo = 21 pI 



Jp2 
d'Eoip) 



dp^ 



P=PO 



P=PO 



sat,2 



K.. 



asy 



^asy „ ^5 

Ksym ~ 6L. 



(13) 
(14) 

(15) 
(16) 



The Eq(p) = E(p, 6 = 0) is the energy per nucleon for symmetric nuclear matter. The is 
the incompressibility coefficient of the symmetric nuclear matter at the saturation density which 
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together with ^sat,2 can yield the value of incompressibility coefficient for asymmetric nuclear 
matter [|l8n . The constraints on the values of ^sym, h ^asy llHSD, IlZLIsI, 
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m, K,,^2 m are, 



£sym = 30 + 5MeV 



L = 88 + 25MeV 



K^,y = -500 + 50MeV 



Ko = 240 + 20MeV 
Ks^t2 = -370+ 120MeV. 



(17) 



IV. ASYMMETRIC NUCLEAR MATTER AND NEUTRON-SKIN 



We study the properties of asymmetric nuclear matter for 26 different parameterizations of the 
ERMF model. Each of these parameterizations are obtained by fitting the experimental data on 
the binding energies and charge radii for few closed shell nuclei. Twenty one of these parameter 
sets were obtained in our earlier work ll25n which correspond to different combinations of neutron- 
skin thickness AR in ^''^Pb nucleus and the co -meson self-coupling strength = ^gl^ (Eq. H]). 
Hereafter, we name these 21 parameter sets as BSRl - BSR21 according to the values of ^ and 
AR. The parameter sets BSRl, BSR2, ...,BSR7 correspond to AT? = 0.16,0.18, ....,0.28 fm with 
^ = 0. Similarly, the parameter sets BSR8, BSR9, BSR14 and BSR15, BSR16, ...,BSR21 
correspond to AR = 0.16,0.18, ...,0.28 fm, but, with ^ = 0.03 and 0.06, respectively. The set of 
experimental data for the binding energies and charge radii used to generate the parameter sets 
BSRl - BSR21 are exactly the same. It can be seen from Ref. [125Q that the rms error on the total 
binding energy and the charge radius for all our 21 parameter sets are more or less the same. The 
rms errors for the total binding energy are 1.5 - 1.8 MeV and that for the charge radii lie within 
the 0.025 - 0.04 fm. Other five parameter sets considered are, FSUGold im, FSUGZ03 [27], G2 
urn, TMl ilJ and TMl* U. 

The density dependence of the symmetry energy E^y^ip) plays central role in understanding the 
behaviour of the asymmetric nuclear matter. The E^y^ip) at subnormal densities can be expressed 
in terms of the £'sym(po), slope L and the curvature ^sym or ^asy as given by Eq. (fT^ . We mainly 
focus on the values of ^symOoo), L and ^asy as obtained for different parameterizations of the ERMF 
model. In addition, we also calculate the curvature parameter ^sat,2 (Eq. [T5l) which together with 
^0 yields the values for the incompressibility coefficient for the asymmetric nuclear matter lll8ll . 
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In Table U we list the values of the AR, B/A, po, ^o, £'sym, L, K^ym, K^sy and ^sat,2 for different 
parameterizations of ERMF model. For the comparison, in the last row of TableHl we also provide 
the results obtained for most commonly used NL3 parameterization [|3O0 of the conventional RMF 
model which includes the non-linear terms only for the cr meson. For better insight, in Fig. [H 
the values of ^symOoo) and L are plotted against AR and ^ for the parameter sets BSRl - BSR7 
(squares), BSR8 - BSR14 (circles) and BSR15 - BSR21 (triangles). Similar plots for Ko, K^sy 
and ^"531,2 are displayed in Fig. |2l The results for the other parameter sets FSUGold, FSUGZ03, 
G2, TMl and TMl* are also depicted in Figs. [T]and[2l We first focus on the results obtained 
for our 21 parameter sets BSRl - BSR21. It can be seen from Figs. [U and [2] that the values of 
the isovector quantities E'symOoo)^ L and ^asy depends mainly on the value of AR. The isoscalar 
quantity Kq is more sensitive to the choice of ^. The values of K^ata^ however, depends on both 
the ^ and AR. The constraints (Eq. [TTl) on the values of ^symOoo) and L are satisfied by our 
parameter sets for which AR = 0.18 - 0.24 fm. The constraint on the values of Kq is satisfied by 
all the 21 parameter sets BSRl - BSR21. The constraints on the value of ^asy are satisfied by our 
parameterizations with AR = 0.22-0.24 fm. On the other hand, except for three parameter sets, all 
of our other parameterizations satisfy the constraint on the value of ^sat,2- These three parameter 
sets are BSR5, BSR6 and BSR7 which correspond to ^ = with AR = 0.24, 0.26 and 0.28 fm, 
respectively. Other parameter sets FSUGold, FSUGZ03, G2, TMl and TMl* are not compatible 
simultaneously with the various constraints summarized in Eq. (fTTl) . In short, it appears that only 
five parameter sets BSR4, BSRl 1, BSRl 2, BSRl 8 and BSRl 9 with AR = 0.22 - 0.24 fm obey the 
empirical constraints of Eq. (flTl) very well. Our parameter sets BSR3, BSRIO and BSR17 with 
AR = 0.20 fm also satisfy all the constraints of Eq. (flTl) except for the ^asy ■ The values of ^asy for 
these parameter sets are only marginally away from the ones given by Eq. (fTTl) . So, on the basis 
of the these constraints, the ERMF model predicts the value of neutron-skin thickness in the ^"^Pb 
nucleus to be ~ 0.20 - 0.24 fm. Similar investigations using the Skyrme Hartree-Fock models 
predicted AR = 0. 18 -0.26 fm. Thus, the empirical constraints extracted from the isospin diff'usion 
data predict more or less model independent values for the neutron-skin thickness. In Ref. [1120, 
various parameter sets of the ERMF model considered are the FSUGZOO, FSUGZ03, FSUGZ06 
and G2. For these cases, AR is either ~ 0.19 or ~ 0.26 fm (see also Table U). These values of 
AR are seem to be either little smaller or quite larger. Consequently, none of the parameter sets 
considered in Ref. [|l2ll are consistent with all the constraints of Eq. (fTTl) . 

The symmetry energy and its density dependence in the ERMF model is mainly governed by 
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values of the coupling strengths gp, rip, rjip and ri2p (Eqs. [21 and (6]). The strengths rjp, rjip and ri2p 
determines the contributions of the mixed interaction terms which account for the coupling of the 
isoscalar cr and oj mesons to the isovector p mesons. It may be emphasized, in the conventional 
RMF models, the contributions of these mixed interaction terms are ignored (i.e., T]p= riip= T]2p = 
0). We have used our parameter sets BSRl - BSR21 to look into the variations of the gp, rip, rjip 
and ri2p with AR. As an illustration, in Fig. [31 we plot the values of the gp/4-n, rjp, rjip and ri2p 
for the parameter sets BSR8 - BSR14 which correspond to different values of AR with ^ = 0.03. 
Scenario for the parameter sets BSRl- BSR7 and BSR15- BSR21 (not shown here) is analogous 
to that of BSR8 - BSR14. We can see from Fig. [3l there is an overall decrease in the values of 
these coupling strengths with increase in AR. It is interesting to note that all the coupling strengths 
are nearly unity for AR ~ 0.22 fm. Coincidently, for AR ~ 0.22 fm, the empirical constraints on 
the behaviour of the symmetry energy as extracted from the isospin diffusion data are also satisfied 
very well. The strengths rjp, rjip and ri2p tend to vanish for higher values of AR. In other words, 
if the contributions of the mixed interaction terms in question are ignored, as in the case of the 
conventional RMF model, the parameters obtained by fit to the experimental data on the binding 
energy and charge radius would give rise to AR > 0.26 fm. The values of AR ~ 0.22 fm, as 
favoured by the isospin diffusion data, can be achieved within the conventional RMF model only 
at the expense of the quality of fit to the bulk properties of the finite nuclei. 

V. ROLE OF MIXED INTERACTIONS 

We would like to investigate the role of various mixed interaction terms which are crucial 
in determining the density dependence of the symmetry energy. In particular, we investigate the 
effects of the terms which account for the coupling of the isoscalar cr and co mesons to the isovector 
p mesons. The coupling constants for these terms are rjp, rjip or rj2p (Eq. [6]). The term with coupling 
constant rjp is of the cubic order in the meson fields. Whereas, the terms with rjip and ri2p are of 
quartic order in the meson fields. Our objective is to delineate the effects of these cubic and quartic 
order mixed interaction terms. For this purpose, we generate two different families of interactions. 
For the first family of interactions Fl, we put rjip = r]2p = and fit the remaining coupling constants 
using appropriate set of experimental data for the bulk properties of finite nuclei. But, the coupling 
constants rjip and rj2p are also included in the fit for the second family F2. Thus, the behaviour of 
the symmetry energy for the Fl family is governed by the coupling constants gp and rjp only. In 
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case of the F2 family, the behaviour of the symmetry energy depends additionally on rjip and rjjp. 
For both the families we obtain the parameter sets corresponding to different values of AR in the 
range of 0.18 - 0.26 fm with fixed ^ = 0.03. The procedure for calibrating the parameters of the 
model as well as the set of experimental data for the binding energies and charge radii used here 
are exactly the same as in Ref. [|27n . In addition, the parameters are also subjected to the empirical 
constraints ofEq. (flTl) . 

In Fig. |4]we plot our results for the rms errors on the total binding energies and charge radii. 
We see that quality of the fits to the total binding energies and charge radii are more or less the 
same for both the Fl and F2 families of interactions, except for AR < 0.2 fm. Therefore, it appears 
that the quartic order mixed interaction terms with coupling constants rjip and ri2p are redundant. 
More precisely, one might say that the values of rjip and ri2p can not be appropriately determined 
by the bulk properties of the finite nuclei. These quartic order terms might play important role in 
fixing the behaviour of the symmetry energy at supranormal densities which is largely unknown at 
present. It can be easily concluded from Fig. |4]that the most preferred value for the neutron-skin 
thickness in the ^°^Pb nucleus within the ERMF model is AR ~ 0.22. In Table UIl we give the pa- 
rameter sets for the Fl family of interactions obtained for different values of AR. These parameter 
sets are named as BKA20, BKA22 and BKA24 which correspond to AR = 0.20, 0.22 and 0.24 
fm, respectively. In Table HIH we present the results for the various quantities associated with the 
symmetric and asymmetric nuclear matter calculated at the saturation density using the parameter 
sets BKA20, BKA22 and BKA24. For the sake of completeness, we have repeated our calculation 
for AR = 0.22 fm with rjp = rjip = rj2p = 0. In this case, the rms errors 6B = 4.8 MeV and 
Srch = 0.05 fm are significantly higher compared to the ones obtained for the Fl and F2 families. 
Thus, the contributions of the mixed interaction terms seem indispensable in order to satisfy si- 
multaneously the empirical constraints on the density dependence of the symmetry energy as well 
as the experimental data on the bulk properties of the finite nuclei. 



VI. CONCLUSIONS 



The density dependence of symmetry energy and the incompressibility coefficient for the asym- 
metric nuclear matter are studied using 26 different parameterizations of the ERMF model. The 
model includes the contributions from self and mixed interaction terms for cr, co and p mesons 
upto the quartic order. Each of the parameterizations considered are compatible with the bulk 
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properties of the finite nuclei. Furthermore, these parameter sets yield the neutron-skin thickness 
in 208p|5 nucleus which vary over a wide range from 0.16 - 0.28 fm. The behaviour of symmetry 
energy at subnormal densities for several parameterizations of the ERMF model corresponding to 
AR ~ 0.20 - 0.24 fm are found to be consistent with the empirical constraints on them as extracted 
from the analyses of the isospin difl'usion data. The ERMF model prediction for AR ~ 0.20 - 0.24 
fm is in reasonable agreement with the ones obtained in the similar way, but, for the Skyrme 
Hartree-Fock model We have investigated the role of the cubic and quartic order mixed in- 
teraction terms which are crucial for the density dependence of the symmetry energy. It is utmost 
important to include the contributions at least from the cubic order term to incorporate the empiri- 
cal constraints on the density dependence of the symmetry energy without affecting the quality of 
the fit to the bulk properties of the finite nuclei. The mixed interaction terms of the quartic order 
might be important to obtain the appropriate behaviour of the symmetry energy at supranormal 
densities which is largely unknown. 
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TABLE I: Some bulk properties of the nuclear matter at the saturation density (po): binding energy per nu- 
cleon {B/A), incompressibility coefficient for symmetric nuclear matter (Kq), symmetry energy (^symOoo)), 
linear density dependence of the symmetry energy (L) and various quantities (A'sym), (^asy) and (^sat2) as 
given by Eqs. (I13H15I ). The values for the neutron-skin thickness AR for the ^°^Pb nucleus are also listed. 
The superscript 'c?' on several parameter sets indicate that they satisfy the constraints of Eq. ([TT] ) very well. 



Force 




B/A 


PO 


^0 


Eavm (Pn) 

tiyy 111 vj / 


L 










(fm) 


(MeV) 




(MeV) 


(MeV) 


(MeV) 


(MeV) 


(MeV) 


(MeV) 


BSRl 


0.16 


16.0 


0. 


148 


240 


31.0 


60 


13 


-345 


-335 


BSR2 


0.18 


16.0 


0. 


149 


240 


31.4 


62 


-4 


-376 


-363 


BSR3 


0.20 


16.1 


0. 


150 


231 


32.6 


71 


-8 


-431 


-395 


BSR4" 


0.22 


16.1 


0. 


150 


239 


33.0 


73 


-20 


-460 


-460 


BSR5 


0.24 


16.1 


0. 


151 


236 


34.3 


83 


-13 


-513 


-509 


BSR6 


0.26 


16.1 


0. 


149 


236 


35.4 


86 


-48 


-562 


-557 


BSR7 


0.28 


16.2 


0. 


149 


232 


37.0 


99 


-15 


-608 


-598 


BSR8 


0.16 


16.0 


0. 


147 


231 


31.0 


60 


-1 


-363 


-286 


BSR9 


0.18 


16.1 


0. 


147 


233 


31.6 


64 


-12 


-396 


-313 


BSRIO 


0.20 


16.1 


0. 


147 


228 


32.6 


71 


-17 


-442 


-361 


BSRll"" 


0.22 


16.1 


0. 


147 


227 


33.6 


79 


-25 


-497 


-387 


BSR12" 


0.24 


16.1 


0. 


147 


232 


33.8 


78 


-44 


-511 


-412 


BSR13 


0.26 


16.1 


0. 


147 


229 


35.6 


91 


-40 


-585 


-466 


BSR14 


0.28 


16.2 


0. 


147 


236 


36.1 


94 


-41 


-601 


-474 


BSR15 


0.16 


16.0 


0. 


146 


227 


30.9 


62 


-22 


-393 


-252 


BSR16 


0.18 


16.1 


0. 


146 


225 


31.2 


62 


-25 


-399 


-258 


BSR17 


0.20 


16.1 


0. 


146 


222 


31.9 


68 


-32 


-437 


-287 


BSRIS'' 


0.22 


16.1 


0. 


146 


221 


32.6 


73 


-42 


-478 


-317 


BSR19" 


0.24 


16.1 


0. 


147 


221 


33.6 


79 


-50 


-526 


-350 


BSR20 


0.26 


16.1 


0. 


146 


223 


34.3 


88 


-39 


-565 


-365 


BSR21 


0.28 


16.1 


0. 


145 


220 


35.7 


93 


-45 


-600 


-402 


FSUGold 


0.21 


16.3 


0. 


148 


229 


32.5 


60 


-52 


-412 


-276 


FSUGZ03 


0.19 


16.1 


0. 


147 


233 


31.6 


64 


-11 


-396 


-312 


G2 


0.26 


-16.1 


0. 


153 


215 


36.4 


100 


-7 


-611 


-404 


TMl 


0.27 


16.3 


0. 


145 


281 13 


36.8 


111 


34 


-632 


-518 


TMl* 


0.27 


16.3 


0. 


145 


281 


37 


102 


-14 


-625 


-429 



TABLE II: Parameter sets for the Fl family of interactions obtained for different values of neutron-skin 
thickness in the ^osp^^ nucleus. These parameter sets are named as BKA20, BKA22 and BKA24 which 
correspond to the neutron-skin thickness 0.20, 0.22 and 0.24 fm, respectively. The couplings T]ip and r]2p 
are taken to be zero. The masses for u) and p mesons are ^ = 782 MeV and ^ = 110 MeV with nucleon 
mass M = 939 MeV. 



Force 


An 


An 


Sp 
An 




^0 


1\ 


m 


1p 


M 


BKA20 


0.8042 


1.0102 


0.9812 


1.1523 1.9892 


4.8344 


0.0005 


0.0657 


3.6164 


0.5430 



BKA22 0.8462 1.1089 1.0302 1.5500 2.13451 5.8253 0.1555 0.0697 3.9294 0.5302 
BKA24 0.8593 1.1463 0.9381 1.7719 3.1064 6.2247 0.2700 0.1159 2.4900 0.5248 



TABLE III: Same as TableU but, for the forces BKA20, BKA22 and BKA24. 



Force 


^R 


BjA 


Po 


^0 


^sym 


L 


^sym 


^asy 






(fm) 


(MeV) 


(fm-') 


(MeV) 


(MeV) 


(MeV) 


(MeV) 


(MeV) 


(MeV) 


BKA20 


0.20 


16.1 


0.146 


240 


32.3 


76 


-15 


-469 


-320 


BKA22 


0.22 


16.1 


0.148 


227 


33.3 


79 


-9 


-483 


-382 


BKA24 


0.24 


16.1 


0.148 


228 


34.3 


85 


-15 


-525 


-420 
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0.16 



0.22 0.24 
AR (fm) 



0.26 



FIG. 1: (colour online) The symmetry energy f'symO'o) and its slope L plotted against the neutron-skin 
thickness Mi in the ^°^Pb nucleus for 26 different parameterizations of the ERMF model (see also Table IJl. 
The open squares, circles and triangles represent the results for the parameter sets BSRl - BSR7, BSR8 - 
BSR14 and BSR15 - BSR21, respectively. The symbol plus with different colours depict the results for the 
FSUGold (green), FSUGZ03 (black), G2 (purple), TMl (maroon) and TMl* (orange). The dashed lines 
represent the constraint on £'sym(po) and L as given in Eq. ([17] ). 
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FIG. 2: (colour online) Same as Fig. \T\ but, for Kq, A'asy and A'sat,2 



16 



12 



-a 

CO 

• 1— ( 

o 

u 



8 - 



4 - 







1 


1 




1 


1 




1 














H Tl 

■ 'p 














♦ ^ip - 




A 










A \p 


— 


♦ 










— 


- 












- 


— ■ 


■ 










— 






■ 















♦ 


■ 






_ 


- 






A 




■A. 


- 

♦ 


_ • 


• 


• 


• 








1 


1 




1 


1 


T 


f 


0.16 


0.18 


0.2 


0.22 


0.24 


0.26 


0.28 



AR (fm) 



FIG. 3: (colour online) The coupling strengths gpl4n, rjp, r]\p and T]2p for the parameter sets BSR8 - BSR14 
plotted against A/?. 
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FIG. 4: (colour online) Plots for the rms errors on the total binding energy 5B (upper panel) and the charge 
radius 5rch (lower panel) for the Fl and F2 families of the interactions. 
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